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ABSTRACT 


A rectangular cross section channel with both straight 
and curved test sections was used to investigate the effect 
of Taylor-Goertler vortices on the heat transfer from the 
curved section by comparing the results with those obtained 
in the straight section. The flow was heated through the 
Joulean heating of one of the walls which was made of 
Temsheet. Liquid crystals were applied on the surfaces of 
both test sections. Experiments were run for several values 
of the volumetric flow rate corresponding to Reynolds num- 
bers between 258 and 920. The average Nusselt number was 
calculated for each experiment and plotted as function of 
the Reynolds and Dean numbers. Three distinct regions of 
the average Nusselt number were observed in the curved 
section when compared with the straight section. The presence 
of Taylor-Goertler vortices in the curved section was veri- 
MmUeGmErone ele Gistribution Of the liguid crystal color 
bands. It was concluded that the overall effect of the 
vortices whenever they occurred was to improve the heat 


transfer process from the curved section. 
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I. INTRODUCTION 


Sete BOR-GORRTLER VORTICES 

Several studies have shown that a fully developed 
laminar flow of a viscous fluid along a concave wall does 
not remain two-dimensional. Instead, the flow forms a 
system of counter-rotating vortices - the Taylor-Goertler 
vortices - whose axes are in the direction of the main flow 
and with secondary velocities in both perpendicular direc- 
tions as shown in Figure l. These vortices are caused by 
the variation in centrifugal forces on the fluid particles 
at different locations in the flow. 

In a curved channel such as used in the present inves- 
tigation the velocity profile of the fully developed laminar 
flow is nearly parabolic. The fluid near the center of the 
channel has a relatively higher velocity and is subjected 
to greater centrifugal forces than the fluid moving near 
the walls. Thus, the tendency is for the fluid in the cen- 
ter of the channel to move outward toward the concave wall. 
The fluid near the concave wall is unable to resist this 
action and must move in the spanwise direction and then 
radially inward replacing the central fluid. Once in the 
center of the channel, this fluid is subjected to the higher 
mew Veloctty and the rotation is continued. This cyclic 


motion forms the counter-rotating Taylor-Goertler vortices. 
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Figure l. Sketch of Taylor-Goertler 
Vortices 





Although the Taylor-Goertler vortices are laminar 
Voretwces ,eteney May affect the transition from laminar to 
turbulent flow and are therefore important to understanding 
this phenomenon. 

The influence of the Taylor-Goertler vortices may also 
explain the increase in the observed heat transfer rates 
from curved walls. 

The Taylor-Goertler vortices have similarities to other 
vortex flows such as the longitudinal vortex rolls developed 
in the laminar forced convection heating of fluid layers 
between parallel plates. The striations seen at stagnation 
points on bluff bodies and the cross hatching observed in 
reentry vehicles have also been explained, at least in 
part, by the presence of Taylor-Goertler vortices. 

There are many possible applications that could result 
from investigation and understanding of the Taylor-Goertler 
vortices. A few of such applications include improved 
cooling of turbine blades and other external surfaces, 
reduction of pressure losses in bends and a better under- 
Staneang GEL the transition from laminar to turbulent flow. 
Heat exchangers could be designed to take advantage of the 
low pressure drop associated with laminar flow and the 
improved heat transfer characteristics of the flow in the 


presence of the vortices. 
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pee ere srl orORyY 

Lord Rayleigh [Ref. 1] in 1916 was the first to consider 
the instability of an inviscid fluid rotating symmetrically 
about an axis. In his study he determined the stability 
Gpueerron fOr ene Inviscid fluid to be that the product 
of the local circumferential velocity and the local radius 
of curvature either increases or at least remains constant 
as the radius increases. 

In 1923 G.I. Taylor [Ref. 2] in an extensive analytic 
and experimental study extended the work of Rayleigh to 
viscous fluids. In his investigation for the flow between 
two rotating cylinders, he found that when the inner cylin- 
der rotates and the outer cylinder is stationary the motion 
of the fluid becomes unstable when the value of the Taylor 


number defined as: 


on d 
Ta = Re ar 
exceeds 41. In the equation above, dis the spacing between 


the cylinders, assumed small when compared with Ri, the 
radius of inner cylinder and Re is the Reynolds number based 
on the circumferential velocity of the inner cylinder and 
the length d. For values of the Taylor number greater than 
the critical value a secondary motion sets in and the 


Taylor vortices occur. 
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A Similar type of instability is observed when a 

viscous fluid flows in a curved channel due to a pressure 
gradient acting around the channel. This problem was first 
considered analytically by W.R. Dean [Ref. 3] for a channel 
formed by two concentric cylinders where the spacing between 
the cylinders was small when compared with the radius of 

the inner cylinder. Dean concluded that the instability 
will first arise and secondary flow vortices similar to the 


Taylor vortices will form when the Dean number defined as 


Z d 
De = Re Ri 
1s greater than 36. In this case the Reynolds number was 


based on the mean velocity of the unperturbed flow and the 
length d. 

The analytical work of Dean was later verified by W.H. 
Reid [Ref. 4] using a simplified method. H. Goertler [Ref. 
5] in 1940 studied the stability of laminar boundary layer 
Peerates On slightly curved walls relative to small distur- 
bances. Goertler found that these disturbances were similar 
to those investigated by G.I. Taylor which led to the 
Peceteiltey In the fFOrm Of vortices. From numerical calcu- 
lations he concluded that amplified disturbances were pro- 
duced only on the concave walls. The approximate results 
obtained by Goertler were verified in 1955 with an exact 


solution developed by G. Hammerlin as reported by H. Schlichting 


7, 
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(-Cromcumomemmanre recently by A.M.O. Smith [Ref. 7] through 
an extensive numerical analysis. 

ies oo, Mae reith (Ref. 8] investigated the influence 
of curvature on heat transfer for fully developed turbulent 
flows and reported that the heat transfer from a heated 
concave wall was considerably higher than from a convex 
wall of the same curvature and under the same conditions 
Gren VOW . 

L. Persen [Ref. 9] in 1965 was able to relate the increase 
in the heat transfer rate from curved walls with the presence 
of the Taylor-Goertler vortices. Persen considered the 
Special cases of a very high and very low Prandtl number and 
determined that the overall effect of the vortices was to 
increase the rate of heat transfer through the boundary 
layer. 

ieee) Pee Meeormack,et al. [Ref. 10] reported the first 
experimental work on the effect of the Taylor-Goertler vor- 
tices on the heat transfer. It was pointed out in their 
Study that in order to explain theoretically the effect of 
the vortices, it was necessary to retain the non-linear 
terms in the flow equations. This was considered by R. 
Kahawita and R. Meroney [Ref. 11] who concluded that the 
high-order terms and the normal velocity components of the 
main flow become increasingly important at small wave 


numbers. 
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Vowviommmanma y. Uchida [Ref. 12] in 1966, considered 
the fully developed forced convective heat transfer between 
horizontal flat plates. They observed that when the tem- 
perature difference between the plates increased above a 
Gritical value, there was a formation of longitudinal vortex 
rolls, similar to the Taylor-Goertler vortices, with the 
axes parallel to the flow direction. M. Akiyama, et al., 
fet. Horne ecoOntirmed the onset of these vortices at 
a critical Rayleigh number of 1708 for the case of heating 
from below. 

The published literature concerning the flow and heat 
transfer in curved rectangular channels has been rather 
limited as reported by R. Shah and A. London [Ref. 14]. 

In 1970, K. Cheng and M. Akiyama [Ref. 15] developed a 
numerical solution for the forced convection heat transfer 
in a curved rectangular channel but only for small aspect 
ratios (nearly square cross section). Analytic and experi- 
mental results for a fully developed constant wall heat 
flux square cross section channel were obtained by Y. Mori, 
emcee. 16). kK. Cheng, et al., [Ref. 17] in 1974 
solved the Graetz problem in a curved square channel using 


a numerical method. 


ro 





ie INTENT OF@EHE STUDY 


The purpose of this study was to investigate the 
effects of the Taylor-~-Goertler vortices on the heat trans- 
fer in a curved rectangular channel and to compare with the 
heat transfer in a straight channel of Similar cross section. 
It was expected that the effect of the vortices would result 
from the secondary velocity components transporting heated 
fluid from the curved wall toward the center of the channel 
and carrying cooler fluid from the center to the curved 
wall. It was also expected that the heat transfer process 
from the curved wall would show improvement as a consequence 
of the presence of the vortices. 

For the purpose of this study, a rectangular cross 
section channel with both straight and curved test sections 
waS considered. The results obtained in the curved section 
were compared with those obtained in the straight section 
to determine the effect of the Taylor-Goertler vortices on 
the heat transfer. 

An analytical solution for the case of the straight 
Beee Section was given by P. McCuen, et al., [Ref. 18] in 
their 1962 study. The results obtained in the straight 
section were compared with the analytical results of McCuen. 

For the curved section, the results obtained in the 


present investigation were compared with the numerical 
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solution developed by K. Cheng and M. Akiyama [Ref. 15] 
for the heat transfer in rectangular curved channels with 


small aspect ratios. 
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IflIfl. EXPERIMENTAL WORK 


A. DESCRIPTION OF THE APPARATUS 

A rectangular cross section channel built as described 
in Refs. 19 and 20 was modified to include two heated 
sections in order to achieve the siiceeere Of “thas Peay 
The channel was made of two 0.635 cm thick sheets of Plexi- 
glas separated by 0.635 cm thick spacers which also formed 
the sides of the channel. 

The channel was composed of a straight portion 122.0 
cm long, a curved section which completed a 180 degrees 
turn and a final short straight section as shown in Figure 
2. The radius of curvature of the interior concave wall 
Spee wecumved pPOortlon was 30.5 cm. The cross section of the 
channel was 0.635 cm high and 25.4 cm wide with an aspect 
ratio of 40 and a cross sectional area of 16.129 square 
centimeters. The cross section of the channel at both test 
sections 1s shown in Figure 3. 

The working fluid was air which entered the channel 
through an aluminum entrance nozzle attached to the 
beginning of the straight portion and which was covered 
with cheesecloth during the experimental work. At the end 
of the channel, the flow was conducted through an aluminum 
exhaust nozzle connected by flexible tubing to a Fisher and 
Porter Company variable area flow meter, model number 


HMOs o6D. this rotometer had a 100% full scale flow rate 
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foo amen ore meters of air per minute (11.1 standard 
cubic feet of air per minute). 

The flow of air was drawn through the channel, the 
Flexible tubing and the rotometer by two electrically driven 
Cadillac centrifugal blowers, model G-12. The blowers 
were connected in parallel, and their speeds were controlled 
by varying the motor voltage with a General Radio Company 
Variac Autotransformer type W10MT3. The voltage was 
regulated by a Sorenson A.C. voltage regulator, model 
ROD 0. 

To obtain the experimental heat transfer data, two 
test sections were constructed: a straight test section 
eee em Ong in the straight portion of the channel after 
the hydrodynamic entrance region and a curved test section 
aaeoecm Long sSsubtending an arc of 53.1° in the curved 
portion of the channel. The area of the straight test 
section was 741.93 square centimeters and that of the curved 
test section was 717.74 square centimeters. 

In each test section the outer wall was removed and 
modified as shown in Figures 4 and 5. A piece of Temsheet - 
a carbon impregnated porous paper with the property of 
uniform electrical resistivity - was glued to the entire 
interior surface of the outer wall in each section. The 
flow of air was heated through the Joulean heating of the 
Temsheet. As the electrical resistance of the Temsheet was 


not constant but slowly variable with the temperature, a 


Ps) 
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Plexiglas \\ 
Temsheet \\3 


Thermocouples \| 





on the 
Temsheet 
Thermocouple 
suspended 
in the flow 
Figure 5. Internal Arrangement of the 


Curved Test Section 
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precision resistor with an electrical resistance of 2.013 
ohms was inserted into the heating circuit in series with 
the Temsheet to calculate the instantaneous power supplied. 
The electrical current to heat the Temsheet was supplied 
from a LAMBDA Regulated Power Supply model LK345A FM. 

The variables to be measured were: the flow inlet 
temperature at the entrance of the channel (Toa) the flow 
outlet temperature after each test section (T ), the wall 
), 


), 


O)UIEE 


temperature of the heated plate in each test section (Tro 


the flow exit temperature at the end of the channel (Toyter 


the temperature of the insulation for each test section 


( ) 


Try) and the voltages across the precision resistor (V 


and the heater (Vis) > 


PR 


Thermocouples were used to measure the temperatures. 
Cholesteric liquid crystals were used also to measure quali- 
tatively the surface temperature distribution on the Temsheet. 
bigmad crystals exhibit dramatic changes in color for small 
changes in temperature. (An excellent reference on the use 
of liquid crystals is the article by Cooper, et al. [Ref. 
21]). The liquid crystals would provide a reasonable 
Gualitative comparison between the heat transfer process 
on both test sections through the observation of the dis- 
PmEaeion of the color bands in each case. The occurrence 
of color stripes on the curved wall, corresponding to areas 
of relatively high and low temperature as a result of the 
presence of the Taylor-Goertler vortices was expected to 


be observed. 
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To select the appropriate liquid crystals to be applied 
on the surface of the Temsheet, the values of the flow inlet 
temperature, the Reynolds number, the Prandtl number and 
the Nusselt number were assumed and an approximate energy 


balance was established: 





Q = m Cp (Toyp- Try) = A Ap, AT 
(ORe 
AT _ Re Pr , Ac) 
Tour ~ Trin ae “Seip 
where 
T +7 
mas ae, . IN OUT 


For a desired difference between the outlet and the 
inlet temperatures of fifteen degrees Celsius, a wall 
temperature of fifty-five degrees Celsius was calculated. 

The required power considering the possible losses was 
estimated between 45 and 75 watts. 

Based on the calculated wall temperature the NCR Micro- 
encapsulated Liquid Crystal R-53 (onset of red at fifty 
three degrees Celsius) was selected. After some preliminary 
experiments, it was found that the response curve of the 


liquid crystal R-53 allowed for a large variation of the 


PLS, 






wall temperature within the range of the color blue. It 

was then decided to use a combination of three different 
liquid crystals instead of a single one in order to obtain 
several color bands and a better appreciation of the average 
wall temperature. A new Temsheet test section was constructed 
and covered with a mixture of the NCR Microencapsulated 
Liquid Crystals R-49, R-56 and S~62 prepared in accordance 
with Ref. 21 and applied directly to the surface with a 
Spray gun. About twenty coats were applied to each plate 

to obtain a clearly observable differentiation between the 
different color bands which it was hoped could be recorded 
photographically. 

Thirty seven copper-constantan 30 gauge thermocouples 
were constructed and inserted at several locations in the 
apparatus to measure the different required temperatures. 
All the thermocouples were connected to a Honeywell thermo- 
couple switch mounted on the side of the apparatus. (See 
Figure 6.) In each test section five thermocouples were 
inserted between the first and the second layers of the 
insulation and another five between the second and third 
layers. These thermocouples were connected in parallel 
to read an average value of the temperature of the insula- 


tion (TL ). Twelve thermocouples, individually connected, 


NS 
were suspended in the flow, four in the entrance of the 
channel, four at the end of the straight test section and 


the last four at the end of the curved test section. At 
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each one of these locations the thermocouples were distributed 
Spanwise at intervals of 5.0 cm. These thermocouples were 
used to measure the bulk temperature of the air. In each 
test section two thermocouples were placed in direct con- 
tact with the Temsheet through holes of very small diameter 
drilled on the middle of each plate at a distance of nO . 2 
cm from the leading and trailing edges. The beads of these 
thermocouples were electrically insulated with ENMAR Heat 
Resisting Glyceryl Phthalate. A final thermocouple was 
inserted into the flexible tubing close to the rotometer. 
All the thermocouples were calibrated using a ROSEMOUNT 
Commutating Bridge, model 920A, and a ROSEMOUNT Calibration 
Bath, model 913A. 

The thermal insulation consisted of ARMSTRONG ARMAFLEX 
22 Sheet Insulation - a flexible foamed plastic material. 
Three layers of insulation, each one with a thickness of 
0.635 cm and covering an area slightly larger than each 
test section, were used. The insulation was fixed in place 
with the use of Adhesive Corrosion Resistant tape as shown 
tier eare 7. 

To prevent any leakage of air into the channel, the 
whole apparatus was sealed with DOW CORNING 781 building 
sealant. 

Two aluminum electrodes 0.318 cm thick were inserted 


between the Temsheet and the plexiglas in each test section 
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General Arrangement of Apparatus and Instrumentation 


Figure 7. 





as shown in Figure 3. The plates with the Temsheet, the 
liquid crystals and the aluminum electrodes were assembled 
in the channel with plastic screws for the purpose of 
electrical insulation. The Thermocouple Switch was con- 
nected to a DYMEC Integrating Digital Voltemeter model 
2401B, and to a reference ice bath. A HEWLETT PACKARD 
Digital Recorder, model 562A, provided an automatic printed 
record of the instantaneous data for each position of the 
thermocouple switch. (See Figure 8.) 

The general arrangement of the apparatus and the 


instrumentation is sketched in Figure 2. 


Be EXPERIMENTAL PROCEDURES 

Experiments were run in each test section for several 
values of the volumetric flow rate in the range between 
0.063 and 0.220 cubic meters of air per minute. 

The appropriate electrical power required for each 
value of the volumetric flow rate in order to heat the 
Temsheet up to a temperature at which the liquid crystals 
would change colors was not known. Several preliminary 
experiments were performed to determine those values of 
electrical power and also to obtain an appreciation for 
the general development of the experiments. 

Estimates of the required time to reach steady-state, 
the necessary data to be taken, the frequency of the 


readings, the overall time of each experiment and the 
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interval between experiments were obtained from the 
preliminary experiments. 

The experiment was considered to have reached steady- 
State when the variation in the average wall temperature 
given by the thermocouples in contact with the Temsheet 
was less than one degree Celsius for a period of ten 
minutes. 

The time to attain steady-state was found to be between 
one hour and one and one half hours. It was decided that 
the duration of all experiments would be two hours. Data 
was taken each half hour ene the first hour and every 
ten minutes in the second hour. 

As the area of the curved test section was slightly 
less than the area of the straight test section, the instan- 
taneous power supplied was different. In order to obtain 
Similar values of dissipated power in both test sections 
for each value of the volumetric flow rate, the power had 
to be adjusted for the curved section during the experiment. 

The room temperature was determined with a mercury 
thermometer located in the vicinity of the apparatus. The 
volumetric flow rate was read directly from the rotometer. 
The instantaneous values of Tw, Toxytp! Tons? TN and Tour 
were automatically printed by the digital recorder in 
millivolts. Based on the results of the calibration of 
the thermocouples, the standard thermocouple tables mene 


used to convert millivolts to degrees Celsius. 


36 








The voltages across the precision resistor (Vip) and 
the heater (V,,) Pe seeseecornaed Inevolts, and thesinstantcaneous 


power supplied (Q,,) was then computed using the expression: 


V 


where Rop was the electrical resistance of the precision 
resistor, previously determined as 2.013 ohms. 

Prior to each experiment, an ice reference bath was 
prepared and used for the single reference thermocouple 
connected in series with the thermocouple switch. 

The steady-state distribution of the liquid crystal 
color bands on the surface of both test sections was 
sketched for most of the experiments and photographed for 
the case of the curved test section using a Nikomat 35 mm 
Gamera wath a Nikon 55 mm £f 3.5 Macro-Nikkor lens. 

During the preliminary experiments the bulk temperature 
thermocouples produced anomolous readings. Determining 
the cause required considerable time. An energy balance 
on the straight test section indicated the value of power 
transmitted by the heater to the air flow to be greater 
tnan the power supplied to the heater. 

It was thought at the beginning that leakage of air 
into the channel could be the cause. However it was later 
shown that the problem consisted of erroneous temperature 


readings given by the thermocouples suspended in the flow. 
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These thermocouples had been inserted into the channel 
through Plexiglas plugs which could be removed for replace- 
ment of any defective thermocouple. 

several tests were performed, and it was determined 
that there was an electrostatic interaction between the 
Plexiglas plugs and the beads of the thermocouples. The 
recorded values of the temperature were from two to five 
degrees Celsius higher than the actual values. This dis- 
tortion was more evident for the higher values of the 
volumetric flow rate. The response time of the thermo- 
couples also was much slower than normal. This was con- 
firmed by removing the thermocouples from the apparatus and 
Observing the approach to room temperature. It was found, 
however, that snae “ine thermocouple beads were removed 
from the vicinity of the Plexiglas plugs they would read 
correct values of temperature with a normal response time. 
As the height of the channel was only 0.635 cm, it was 
difficult to move the thermocouple beads away from the 
Plexiglas plugs. 

Considering that the problem probably resulted from 
some electrostatic charge induced in the Plexiglas and 
transmitted to the beads it was decided to try different 
materials for the construction of the plugs. Phenolic 
resin fiberboard and pine wood were used, but the same 


unsatisfactory results were obtained. 
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Finally it was decided to devise a process to insert 
the thermocouples into the channel without the use of the 
plugs and also to bend the thermocouples in order to place 
them parallel to the air flow as shown in Figures 4 and 5. 


When this was accomplished, the correct data was obtained. 
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IV. PRESENTATION OF DATA 


A. ANALYSIS 

Since Temsheet has the property of uniform electrical 
resistivity, the heated wall was considered to be a constant 
heat flux surface. Also because of the large aspect ratio, 
the channel was considered to be approximately parallel plates. 
It was therefore assumed that the experimental configuration 
represented forced convection between parallel plates with 
one surface subjected to constant heat flux and with some 
heat loss at the opposite surface. To analyse this situation 
several quantities were defined and calculated as follows: 

The average heat transfer coefficient between the heated 
wall and the flow of air in each test section was defined 
by the equation 

Caulie = jal oe AT 

In this equation, Qair was the heat convected- (See the 
Figure in Appendix B.) h was the average heat transfer 
coefficient, Apr, was the area of the heated surface and 
AT was the difference between the average wall temperature 
(Two) and the bulk temperature of the fluid (T,). The bulk 
temperature was defined as the arithmetic average of the 
fluid inlet temperature (Tx) and the fluid outlet temperature 


Gy y. In mathematical form AT was expressed by: 


OUT 






IN OUT 


The heat convected (Qair) was calculated with the 


equation 


Qair = th Cp (Loup - Te) 
where Cp was the specific heat of air at constant pressure 
and m was the mass flow rate of air given by the product 
of the volumetric flow rate (q) and the density of the air 
(0). This density was computed assuming perfect gas 
behavior at atmospheric pressure and using the flow exit 


Ve 


temperature Coyit 


The average Nusselt number was then computed from the 


average heat transfer coefficient through the expression: 





where d was the height of the channel and Kair was the 
thermal conductivity of the air at the bulk temperature 


ae). 


B 
The Reynolds number was calculated for all the 


experiments and defined as: 





Ps 
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where 


Ac was the cross sectional area of the channel and u was 
the dynamic viscosity of the air at the flow exit temperature 


(T ). 


j2 > iGah 
For the experiments in the curved test section, the 


Dean number defined by: 


was also calculated. In this equation, d was the height 
of the channel and Ri was the radius of curvature of the 
interior surface of the convex wall. 

The heat lost by conduction through the outer plate 


and the insulation was determined with the expression: 


QO — 
Le (OXI NS/King Apy)? 


where AX Ns was the thickness of each layer of insulation, 


Kons was the thermal conductivity of the insulation and 


AT oS was the difference in the temperatures given by the 
thermocouples located between the first and second layers 


and those located between the second and third layers. 
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In order to verify that Qair as measured by the 
temperature rise in the air was actually the convection 
heat transfer from the Temsheet, an estimate of the radiation 
Meetesthansper in the channel had to be made. To facilitate 
this, one experiment was run in which the temperature of 
the lower plate was measured with two thermocouples. For 
this experiment the radiated heat transfer was computed 
with the expression: 

) 


fe) (Two" - Twl 


R 


Q 
R 


where o is the Stefan-Boltzmann constant and Rp was the 
total resistance to radiation heat transfer between 


the plates given by: 


In this expression a is the emissivity of the interior 
Surface of the outer plate (Temsheet) and Ring is the 
emissivity of the interior surface of the inner plate 
(Plexiglas). 

A sketch of the control volume for the energy balance 
indicating the different components involved and the sample 
calculations, for the case in which the radiated heat 


transfer was also computed, are given in Appendix B. 
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Bee one oULTS 

The data collected from the experiments was evaluated 
through the expressions described in the previous section 
as illustrated in Appendix B, and the results for the major 
variables involved are shown in Table I for the straight 
test section and in Table II for the curved test section. 
A plot of the natural logarithm of the Reynolds number 
versus the natural logarithm of the average Nusselt number 
moegiven in Figure 9 for both test sections. The incertainty 
bands are shown for two representative points. Sample 
calculations for obtaining these incertainty bands are 
Shown in Appendix A. Figure 10 is a plot of the natural 
logarithm of the Dean number versus the natural logarithm 
of the average Nusselt number for the curved test section. 

The results indicated an increase in the heat transfer 
rates with increasing values of Reynolds number for both 
test sections. However, when comparing the plotted values 
of the average Nusselt number for the curved section with 
those of the straight section, three distinct regions were 
detected in the curved section. For the highest values 
of Reynolds number corresponding to the highest values of 
the volumetric flow rate, the Nusselt numbers for the 
curved test section were greater than those of the straight 
test section. For the lowest values of Reynolds number, 
the opposite was noticed with greater heat transfer rates 
occurring in the straight test section. A transition region 


was observed for intermediate values of Reynolds number. 
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Figure 9. Natural Logarithm of Reynolds Number 
versus Natural Logarithm of Average 
Nusselt Number for Both Test Sections 
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Figure 10. Natural Logarithm of Dean Number versus 
Natural Logarithm of Average Nusselt 
Number for Curved Test Section 
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These differences in heat transfer behavior for the 
curved section as compared with the straight section were 
then related to the occurrence of Taylor-Goertler vortices 
through the observation of the distribution of the liquid 
crystal color bands on the surface of the curved section 
for the different values of Reynolds number. It was observed 
that for the smallest values of Reynolds number the differ- 
ent color bands occurred in almost straight lines perpendicu- 
lar to the direction of the flow. The same had been verified 
Mpeweeseraight section independently of the value of Reynolds 
number. It was then observed that as the Reynolds number 
was increased for the curved section, color stripes oriented 
in the streamwise direction, became visible in the region 
close to the trailing edge. These color stripes were ex- 
plained by the presence of the Taylor-Goertler vortices 
and, for the highest values of the Reynolds number, could 
be seen over the entire surface of the curved test section. 
Schematic drawings of the distribution of the color bands 
Sere Surtace of the straight and curved test sections are 
Shew? in Figures ll, 12, 13 and 14. 

The occurrence of color stripes and consequently Taylor- 
Goertler vortices at the same values of the Reynolds number 
for which the average Nusselt numbers were greater in the 
emieved test section led to the conclusion that the effect 
of the vortices was to increase the heat transfer rates 


Enemmthe curved section. 
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The plot of the natural logarithm of the Dean number 
versus the natural logarithm of the Nusselt number for the 
curved section was then analyzed, and straight lines were 
fitted through the results in each region using least 
Squares regression. These are shown in Figure 10. The 
slopes of these lines for the regions of highest and lowest 
values of Dean number were determined as respectively 0.557 
and 0.591 while in the transition region the slope was 
1.57. The limiting values of Dean number for the transition 
region were calculated from the plot as 57.1 and 67.4. 
Three correlations expressing the Nusselt number as a func- 
tion of the Dean number for each one of the regions were 
determined as follows: 

1. Low Values of Dean Number (De < 57.1) 


Nu = «0.292 pe?°?7+ 


Pee iransition Region (57.1 < De < 67.4) 


a OnOOss6 Demuaee 


3. High Values of Dean Number (De > 67.4) 


Siem ses sales! 
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A single correlation for the overall range of Dean 


numbers was also determined as: 


Sa Oe Der eto 


De 





V. DISCUSSION AND CONCLUSIONS 


The measurement of the surface temperature of the 
unheated plate in the straight test section showed that 
this temperature was very close to the bulk temperature 
eeethe fluid. 

It was then assumed that in the absence of significant 
difference in temperature between the inner wall and the air 
flow, the heat convected to the air was exclusively from 
the heated wall. Based on this assumption it was decided 
to compare the results of the present study for the straight 
test section with the analytical solution of the same prob- 
lem developed by McCuen, et al. [Ref. 18]. The hyarautic 
diameter of the channel was calculated as Dh = 1.239 cn, 
and the values of Reynolds number based on the hydraulic 
diameter as used by McCuen also were determined. A dimen- 


sionless axial length coordinate x, given by 


= af. ee 
Dh Re Pr 


where Pr was the Prandtl number, was calculated for 

x = 29.2 cm (length of the straight test section) for all 
the experiments. Based on the values of x it was determined 
that the development of the experiments would occur essen- 
tially in the thermal entrance region, and it was decided 


to integrate the expression of the local Nusselt number 
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for small values of x 


Nu = al 


a re 
Meenogses - = (en 


Qwo 


Seven by’ McCuen, et al., in order to obtain the following 


expression for the average Nusselt number. 


—— 3 2 Qwi,—2/3 
Nu = = —————— ln[{l - ~———=--— (1 + =—)x ] 
4(1 + QW1) > 0.67095 Qwo 
Qwo 


In this expression Qwi was the heat convected to the air 
from the unheated wall, and Qwo was the heat convected 
from the outer Temsheet wall. The average Nusselt numbers 
were then calculated with the expression developed above 
where, in accordance with the previous assumption, the 
value of Qwi was taken as zero and the value of Qwo was 
considered to be the same as Qair. 

The average Nusselt numbers obtained from the experi- 
mental data in the straight section were recalculated as 
function of the hydraulic diameter and were plotted versus 
Reynolds number as shown in Figure 15 along with the 
analytical expression from McCuen. The only significant 
difference observed for the lowest value of Reynolds 
number may be explained by the fact that the value of x 
was equal to 0.064. In this case the Nusselt number 


expression for small values of x is probably not valid. 
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A comparison of the data from the curved section with 
the analytical study by Cheng and Akiyama [Ref. 15] was 
attempted also for the curved test section, but the results 
were less conclusive. Their investigation was done for 
small aspect ratios (nearly square cross section), and the 
results were plotted only for small values of Dean number. 
An extrapolation of these results to greater values of Dean 
number for the case of the largest aspect ratio (5 to 1) is 
shown in Figure 16 simultaneously with the results of the 
present study. 

After these comparisons, several conclusions about the 
heat transfer in rectangular channels were developed from 
the results obtained in the present investigation. The 
heat transfer rates were demonstrated to increase with 
increasing values of Reynolds number in both test sections. 
The values of the average Nusselt number for the curved 
section when compared with the values for the straight 
section seemed to indicate that there are three distinct 
regions of Reynolds or Dean number which can be considered. 
These regions were: an initial region for the lowest 
values of Dean number where the heat transfer rates from 
the straight test section were greater than from the 
curved test section, a final region for the highest values 
of Dean number where the opposite was true and a transition 


region for the intermediate values of Dean number. 
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The observation of the presence of Taylor-Goertler 
vortices at values of Dean number for which the heat 
transfer rates were greater from the curved test section 
led to the conclusion that the overall effect of the 
vortices was to improve the heat transfer process from the 
curved section whenever their presence occurred. The 
critical value of Dean number for the onset of the vortices 
was determined from Figure 10 as 57.1. 

Finally it seemed apparent from the comparison with 
the work of Cheng and Akiyama that for values of aspect 
ratio greater than five the results are independent of 


the aspect ratio. 
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VI. RECOMMENDATIONS 


Several recommendations can be made to improve the 
experimental study of heat transfer in curved as well as 
straight rectangular channels. More experiments should 
be performed in the transition region of the Dean number 
BOreEne Curved test section. This would permit a better 
understanding of the variation of the average Nusselt number 
with the Dean number and a more rigorous estimate of the 
critical value of Dean number above which the Taylor- 
Goertler vortices will form. The effect of the vortices 
on the heat transfer from the curved section could then be 
explained with more detail, and improved correlations for 
the average Nusselt number as a function of the Dean number 
could be developed. 

With respect to the apparatus, it became apparent that 
an improved method of suspending the thermocouples in the 
flow to measure To and Tour should be devised. Considerable 
thought must be given to this problem in order to eliminate 
the interaction between the Plexiglas and the beads of the 
thermocouples and also to make the distance between the 
bea@s and the plates more uniform for all the suspended 
thermocouples. The surface temperature of the inner wall 
should be measured in future experiments. 

Panally; a motion picture study of the onset of the 
Taylor-Goertler vortices in the curved test section, uSing 


the liquid crystals, would be possible and of great interest. 
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APPENDIX A 


ERROR ANALYSIS 


The uncertainties for the major variables in the 
experiments were calculated in accordance with the method 
Peeaned by S. Kline and F. MeGlancoch [Ret 52212. tne 
estimates of the uncertainty in the measured quantities 
were made quite conservatively so that there was considerable 
confidence in the calculated uncertainties. AS an example, 
the calculation of the uncertainty in the Reynolds number 


is given below. The Reynolds number was defined by the 


equation: 





The uncertainty in the flow velocity was determined as 


I= 0.02176 from similar calculations based on estimates 


U 
of the uncertainties in the cross sectional area of the 
channel and in the reading of the volumetric flow rate. 


The uncertainties in the density, the height of the channel 


and the dynamic viscosity were obtained from estimates as 


oS 





do dd L 


= 0.00082, S* = 0.02000 ana = as 0.002230 0The 
uncertainty in the Reynolds number was then calculated 


as 


oe ~ [(0.00082) + (0. 02176) + (0.02)*+ (0.00223) * = 0.02965 


The uncertainty in the Nusselt number was computed 


with the expression: 
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The values of the uncertainties for other variables 


are listed below. 


Quantity Uncertainty 
Ac 0202010 
Ant OR 00265 
Cp 0.00415 
De 0203576 
h 0.06242 
Kair 0.00038 
Nu 0.06555 
g 0.00835 
Qair 0.04050 
Re 0702965 
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Quantity 
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Uncertainty 


0.01640 
0.00453 
0.01060 
0 032107 
0.03940 
0.04743 
0200223 


0.00082 





APPENDIX B 


SAMPLE CALCULATIONS 


A sketch of the control volume for the energy balance 
on the straight test section indicating the major heat 
transfer components involved is shown in Figure 17. A 
Similar diagram could be considered for the curved test 
eeeceten. Sample calculations of the power stippYied (Qp), 
the heat convected to the air (Qair), the heat lost through 


the insulation (Q._), the radiated heat transfer (On), 


Lo 
the average heat transfer coefficient (h), the average 
Nusselt number (Nu) and the Reynolds number (Re) are given 


below for the straight test section. 
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ENERGY BALANCE 


Temsheet 


Figure 17. Energy Balance in Straight Test Section 


SAMPLE CALCULATIONS 


A. DATA 

= ° 
T ROOM oo. G 
& = 0,126 m°/min 
Vor =——sec25 V 
Vay = 35.580 V 
Ror = (25 Muse 

= : = . ay 
ToNS1 1.997 mv 2S Ls 

—_ Z = 38.33 st 
To s2 i542 mv 

= =" 247258 -G 
TN 0.963 mv 

= ; OPO 2 Gy ee 
Ty2 0,962 mv 
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IN4 


EALIT 
Twol 


Two2 


TouT1 
OUT2 
OUT 


OUT4 
‘Twa. 1 


Dye 


Kalir 
Cpair 
ualr 


INS 


axTNS 


e955 my = 24.10 °C 


Boom. =ie2d.34 °C 
SOV Zomvy = ~ 26.97 °C 
-684 mv = 64.93 °C 
woo cum. = 59219 °C 
Scie = 99 39752 °C 
SOW Genie =. 39.287. °C 
Pocveme = 37,97 °C 
mooztmy~yy = 41.21 °C 
oe ee SD . 345° C 
woe = SI 660° 
PO?4Z m* 

.0016 m@ 

-OUGS5 m 


—OzZoo2 W/m °C 


0057 Kd/Kg °C (0.241 BTU/lbm °F) 


5 


.983 x 10 ~ Kg/m-sec 
nella s ip 4 ives W/m °C 
.00635 m 

90 

.75 

0 


1669 x 100° W/m-°K 
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B. 


TEMPERATURE CALCULATIONS 


1. Average Inlet Temperature (Tay) 
T + + 
Meee GIN aN2 “m3 INA 
IN 4 


Bae ene 272410 + 24.0354 
2 a 


= e425 eC 


2. Average Outlet Temperature (Ty ) 


ME 
+ 
, _ *ouri * Tour2 * Tours * Toura 
OUT 4 
B9no2 + 39.87 + 37.97 +. 41021 
4 
=a O64) °C 
3. Temperature Difference in Insulation (AT wg) 
= = = = — ° 
AT SS = TINS TNS? 49.13 Seess 10.80 °C 
4. Average Temperature of Outer Plate (To) 
Tol * Two2 _ 64.93 + 59.19 
it = ie eee = ——————— 
wo 2 2 
aew62.06 °C ol= Soov06e kh 
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>. Average Temperature of Inner Plate (Ts? 


. Cee OT 35394 4 3166 
wi 2 = 2 


Soe oO ete ~Sr 306.50 °K 


6. Bulk Temperature (T,) 


Bo) em eer | 24.25 + 39.64 | : 
3 oe Se =. 8c 


7. Mean Temperature Difference (AT) 


es Oe en =" 62,06e——sl2 905 = 30.11. °C 
wo B 


POWER CALCULATIONS 


1. Power Supplied (0, 


PR ‘H _ (35.580) (3.223) 


2. Heat Lost Through Outer Plate (O, 6) 


oaanG 


QO = 
Lo AXins/Kins “py 


7 ORT U = a a lg a 


(0.00635)/(4.18 x 107°) (0.0742) ) 
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3. Heat Radiated (Q,) 


a. Radiation Resistance (Rp) 





1 - a 1 la aes 
Re = & + + “ 
PL “wo PL wo-wil PL Cwi 
1 : 
oO ak esate | ell een 
PL wo Wl PL 


b. Heat Radiated (O,) 


= a 


a CT =e eee) 
Q = wo wl 
R Re 
-8 4 4 
CS e669 LOns 20742) (3355 bt = S06. 5) 
i444) 
= 2.025 
4. Heat Convected to the Air (Qair) 
a. Volumetric Flow Rate 
b. Density 
So 
Oe ee eens eee 
5“ RE ORD ae 


1.177 Kg/m? (0.073 lbm/ft>) 


fal 





Gs. Heat Convected (O __) 
— Bukig 


tO 
tt 


air Pd CP (Toye 7 Try)? 


(ele On O02 09) ( LO05S7 (2964 24.2 5) 


0.03825 KJ/sec 


Bo, 20 W 


AVERAGE HEAT TRANSFER COEFFICIENT 


Fh Shibsg eo 38.25 
AT (050722) (30-11) 


17.12 W/m* °c 


AVERAGE NUSSELT NUMBER 


—= hd _ (17.12) (0.00635) 


aes (0.02662) Sere 





REYNOLDS NUMBER 
1. Mean Velocity (U) 
q OF00Z209 


AS = (moO LGn = ie 3062Z> m/sec 


CG 
Il 


v2 





2. Reynolds Number (Re) 





ne = 2UA _ (1-177) (1.30625) (0.00635) 
Pea alone 
= 492.33 
Re = 492 


Us 
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